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Abstract 

We study four dimensional large-N SU(iV) Yang-Mills theory coupled to 
adjoint overlap fermions on a single site lattice. Lattice simulations along 
with perturbation theory show that the bare quark mass has to be taken 
to zero as one takes the continuum limit in order to be in the physically 
relevant center-symmetric phase. But, it seems that it is possible to take the 
continuum limit with any renormalized quark mass and still be in the center- 
symmetric physics. We have also conducted a study of the correlations 
between Polyakov loop operators in different directions and obtained the 
range for the Wilson mass parameter that enters the overlap Dirac operator. 
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The large N limit of gauge theories has many intriguing properties. One 
of these is continuum reduction p]. It states that one obtains correct infi- 
nite volume zero temperature results by working on a finite volume lattice 
as long as the center symmetry is intact. In [2], it was proposed, that for 
a Yang-Mills theory with massless adjoint fermions with periodic boundary 
conditions, the volume can be reduced down to a single site as opposed to 
the pure gauge case [3J, where weak coupling analysis shows all the cen- 
ter symmetries to be broken [3]. This has been confirmed both by lattice 
techniques and by perturbation theory [5l[6l[7l[8ll9l [T0l[TTl[T2l[T3] . 

The question we want to address in this paper is what occurs at the large 
N continuum limit when fermions have a mass. The large N continuum 
limit is taken by first extrapolating N — > oo and then b — > oo, where b is the 
inverse 't Hooft coupling, • It has been argued in [13] that for any finite 
mass, a center symmetry unbroken phase exists at sufficiently small volume. 
Lattice studies using Wilson fermions has shown a large range of masses at 
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fixed lattice spacing where the center symmetry remains intact |lL)j . 

In this paper, we address the question of center symmetry both in the 
lattice and in the continuum using massive adjoint overlap fermions |16l 
I17j . We show that the critical bare quark mass [a c , above which the center 
symmetry is broken, is zero at the continuum limit. However, on a lattice 
with a finite lattice spacing, fj, c > 0. Values of masses, which are accessible 
to lattice simulations, depend on how fj, c scales as a function of the lattice 
spacing. 

We study the problem with one Weyl fermion, / = 0.5, both by pertur- 
bation theory and lattice simulations. Using perturbation theory we show 
that center symmetry is broken even when quarks are given an arbitrarily 
small mass. We have performed lattice simulations with different b and N. 
The lattice results confirm with perturbation theory and we find a [i c (b) 
that decreases as b increases. We do not see any evidence of scaling of /i c (b) 
versus b. Our numerical results indicate that we can obtain the continuum 
limit with arbitrary physical mass for the adjoint quarks. 

All details pertaining to the single site lattice model with adjoint overlap 
fermions are described in [9]. To study the continuum limit starting from 
the single site action, we use the weak coupling expansion and write the link 
matrices as 



D 13 



(1) 



and perform an expansion in a^. The 9^ are the eigenvalues of the Polyakov 
loop operator and they have to uniformly distributed in the range [—it, it] 
and uncorrelated in all four directions in order to correctly reproduce infinite 
volume continuum perturbation theory. The leading order result is 
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where the first line of RHS is the contribution from gauge fields [1] and the 
second line is the contribution from / flavors of Dirac fermions [9] . The bare 



quark mass is 



2m w fi 



with /i G [0, 1] and m w is the Wilson mass parameter. 



The gauge action has its minimum, — oo, when all the angles 9^ are equal. 
With one massless Weyl fermion (/ = 0.5) the fermionic part cancels out the 
infinity and renders the action finite. In [9] we used Monte Carlo techniques 
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to find out the actual minimum. Namely, we consider the Hamiltonian 



H = \Y,{<? + P S - ( 3 ) 

For large j3, the Boltzmann measure e~ H is dominated by the minimum. 
Hence, this minimum can be found by performing a HMC update for the 
7r, 9 system. 

To reduce rounding errors in equations of motions, we introduce a regu- 
lator A to the gauge field action 
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In the computations we choose A = 10 -4 , which is much smaller than the 
average difference between angles 2-k/N when N < 200. 

A choice for the order parameters associated with the symmetries 

ism 

If = n, then the Zjv symmetry in that direction is unbroken. 

In Fig. [T] we reproduce the results of [9] with a large N scaling. To better 
observe the center symmetry breaking, the measurements Pi are ordered for 
each configuration s.t. Pi < P2 < P3 < Pa- This indicates that center 
symmetry is probably restored when Wilson mass is in the range 3.0 < 
m w < 10.0. The range of m w with broken symmetry does not depend on N. 

It is possible that center symmetry is broken in a subtle manner in 
the range 3.0 < m w < 10.0. For example, the eigenvalues of the individ- 
ual Polyakov loop operator might be uniformly distributed but they might 
show correlations in different direction^] In order to have a correct sum 
over all momenta as one would have in a infinite lattice, we need to ensure 
that the traces of Polyakov loops vanish and there are no correlations be- 
tween different directions on the single site model. Let us assume that the 
eigenvalues are uniformly distributed and choose ^ff-, j = 1, • • • ,N as the 
TV eigenvalues. Let tTj, j = 1 • • • , N denote a permutation of j = 1 • • • , N. 
We compute the correlated action, S c , with 9^ = and compare it to the 

uncorrelated action, S u , with 9^ = N J where ir^ are different permutations 
for different [i. 

Fig.[2]shows the difference, S u — S c , as a function of Wilson mass m w with 
different N. The value for S u is obtained by averaging over several different 



lr This is the problem with quenched Equchi-Kawai model [15] . 
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Figure 1: Plot of P M as a function of Wilson mass. 




Figure 2: The difference between actions with correlated S c and uncorrelated eigenvalues 
to each direction. 
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Figure 3: Plot of P M as a function of m w with massive quarks. 



random permutations but the fluctuations get smaller as N increases and 
it is sufficient to consider just one random permutation as N — > oo. The 
uncorrelated minimum is preferred when 1 < m w < 5. There is again 
virtually no dependence on N. This combined with the restriction of center 
symmetry restoration gives 

3 < m w < 5 (6) 

as the range for Wilson mass. 

One might wonder why the region of allowed m w does not include zero. 
In a typical free field analysis of overlap fermions |18|. [T9] , one shows that it 
correctly represents a single Dirac flavor in a region around zero momentum 
as long as < m w < 2. Momentum in our case is replaced by (0* — Ofy 
and we want to cover the whole range of allowed momenta. If this does 
not occur, we will not have proper reduction or a correct realization of the 
center symmetric phase. Because the range of allowed momenta (volume of 
the Brillouin zone) in the conventional free field analysis increases as m w 
increases, we see why m w close to zero is not appropriate. Since we do 
not have a concept of doublers on a single site lattice, we do not require 
m w < 2 [9\. Therefore, to reach momenta close to ir and have proper 
sampling of all momenta as per the infinite lattice, we find a range of allowed 
m w than includes m w > 2. One can also understand why m w cannot be 
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Figure 4: Plot of Pi as a function of mass for three different b with N = 15 



arbitrarily large since we would be approaching the limit of naive fermions 
which does not have a center symmetric phase on a single site lattice [9]. 

Once fermions have a non-zero mass, the fermionic contribution to ^ 
is always finite. Then the minimum of the perturbative action is dominated 
by the pure gauge part and occurs when all the eigenvalues are the same. 
The effect of finite N is is demonstrated in Fig. [3] for /i = 0.1. We have only 
plotted the component Pi, since it determines the center symmetry breaking 
point. The symmetry breaking is evident as N — > oo. 

For the actual lattice simulation, we used HMC-algorithm described 
in [9]. All the simulations were performed with / = 0.5, Wilson mass 
m w = 5,0 and they consist of about 100 independent measurements. Ther- 
malization is fast and requires only about ten iterations. Most of the simu- 
lations were performed with N = 15, but to study 1/N effects we did also 
simulations with N = 11 and N = 18. The purpose of the simulations are 
to find out the critical mass \x c for center symmetry breaking as a function 
of N and b. 

In Fig. [4] we have plotted Pi as a function of mass with N = 15 for 



2 This value is slightly high, since it is on the high end of jfi]). This is because the 
argument presented with regard to Fig. [2] was realized after we obtained the numerical 
results presented in this section. 
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Figure 5: Plot of Pi as a function of mass for three different N with 6 = 5 



b = 1,3, and 5. The data shows that /x c (6) does decrease with increasing b 
but the decrease is clearly slower than scaling would dictate. The range of 
center symmetry breaking is between 0.1 and 0.3 for b in the range [1,5]. 
Ignoring wave function renormalization, the dominant part of the scaling 

8tt 2 i, 

dictates that we need to keep iie^ fixed as we take b — > oo in order 
to take the continuum limit at a fixed physical mass. Our data for fJ> c (b), 
therefore, clearly indicates that we can take the continuum limit of a massive 
adjoint fermion coupled to a large N gauge field without any restriction on 
its physical mass. 

To understand the effects of finite N, we performed simulations with 
6 = 5 also at N = 11 and N = 18. The 1/N effects are rather small except 
in the region of the phase transitions as can be seen in Fig. [5j The critical 
value at b = 5 is about 0.1. 
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